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Abstract. Accretion discs are composed of ionized gas in motion around a central object. Sometimes, the 
disc is the source of powerful bipolar jets along its rotation axis. Theoretical models invoke the existence of a 
bipolar magnetic field crossing the disc and require two conditions to produce powerful jets: field lines need 
to be bent enough at the disc surface and the magnetic field needs to be close to equipartition. The work of 
Petrucci et al (2008) on the variability of X-ray binaries supposes that transitions between pure accretion 
phases and accretion-ejection phases are due to some variations of the disc magnetization. This rises the 
problem of the magnetic field dragging in accretion discs. We revisit the method developed by Lubow et al 
(1994) by including momentum and mass conservation equations in a time-dependent ID MHD code. 
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1 Introduction 



First works on magnetic field dragging in accretion discs have focused on the Blandford and Payne (1 19821) neces 



sary c o ndition for cold jet launching, nam ely: tan~^{Br/Bz) > 30. Thus. Ivan BallegooiienI ( 1989 ) . Lubow et al 
(1994), Reves-Ruiz and Stepinskil ( 19961) developed different methods to address this question in a standard ac- 
cretion disc (here after SAD). They all conclude that a large effective Magnetic Prandtl number (here after 
Pm = v^/vm as the ratio of the turbulent v iscosity to the magnet i c diffu sivity) is necessary to reach a sufficient 
magnetic flux accumulation that satisfy the iBlandford and Payne (Il982h criterion. Indeed, in a SAD, bent field 
lines requires a magnetic Reynolds number i?,„ = rUr/vm = RePm ^ T/h (with Re the Reynolds number), 
whereas in a SAD, Re^ 1. This is problematic as it t ends not to be s atisfie d according to 3D MRI simulations 
( Lesur and Longaretti 20091 ). More recent works, by iLovelace et al. I (l2009l) . showed that significant advection 
of magnetic field could be possible even with Pm ~ 1. 



On th e other hand, the existence of jets do require the presence of a large scale vertical field (|McKinnev and Blandford 
20091). It has been shown that this field must be close to equipartition in order to drive powerful jets 
(jFerreira and Pelletieiill995l ). 

In this work, we revisit the standing problem of field advection in accretion discs, by introducing new physical 
input from modern works. 



2 Analytical description 



The full MHD equations could be written in axisymetry in the in the (r, 0, z) cylindrical coordinates. Writting: 
B and b = rB^,, we have: 



a = rA^ with V x A 



t + v.p.^ = o 
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Where 77^ = ^o^-m 

We focus on the magnetic flux transport coupled with mass and momentum conservation. We assume a thin 
keplerian disc with e = h/r constant with radius {h the height of the disc at a given radius). We also assume 
that only the turbulent viscosity plays a role in the angular momentum transfer. 



2.1 Induction equation 

We start with the induction equation, by setting ~ J adz, we rewrite equation 12.31 integrated over the 
thickness of the disc: 

With Uq = Ma/{2TTrY,) with S the surface density, Ma the accretion rate and J^s the surface current density. 
All the question is determining J^g in term of the radial distribution of Approximating the disc as rings of 
toroidal current, an equivalence between ^'(r) and J^g can be found when assuming a potential magnetosphere 



( Lubow et al.lll994f ) 



2.2 Momentum conservation 

By considering a quasi Keplerian motion, one can deduce the accretion rate from equation 

Ma^^^i^^^v^KT^) (2.5) 

iIkt or ^ ' 

2.3 Mass conservation 

We then use the mass conservation fequation l2.1l) integrated over the thickness of the disc: 

9S 1 dMa 



dt 27rr dr 

2.4 Turbulent prescription 

At last, we need to do an assumption on the transports parameters: 



(2.6) 



(2.7) 
(2.8) 



Magnetic fields in accretion discs are known as a potential source of turbulence via the Magneto-Rotational 
Instabilit y (here after MRI). This is e xpected to lead to anomalous transport coefficients Vv and Vm- Recent 
results by lLesur and Longarettil ( 20091 ) show that P,„ ^ 2 and a behaviour of alpha in pP-^ for low (with the 



magnetization: ratio of the magnetic pressure on the thermal pressure). As we know that MRI is suppressed 
around equipartition, we use the prescription: 

a„(M) = 2(^^(1 - /i))i/2 (2.9) 
Where the local disc magnetization is computed as: 

One can study the evolution of the magnetic field and the density in the disc, by solving equations 12. 4[ 12. 6[ 
using [2751 and the equivalence between J^s and '^(r). 
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3 Results 



We study the behaviour of a SAD with Pm = 1 initially thread by a magnetic flux very compressed in the inner 
zone of the disc. This field is chosen t o have a uniform mag netization at t = 0. This initial condition could 
model the high soft state (labelled C in IPetrucci et al.l ( 2008 )). right after the jet quenching of X-ray binaries. 
We use the a„ (/i) prescription (equation 12.91) . The results are shown on figure [T] At the beginning, of the 
simulation, the magnetic field diffuse outward as expected in a SAD. Leading to an increasing of /i in the outer 
zone of the disc, and so a noticeable positive variation of the accretion rate at the outer edge of the disc. We 
can notice that at the end of the simulation, the disc is approaching a stationary state with an accretion rate 
quite uniform at the outer radii of the disc. But as the magnetization is very weak in the inner zone, so is a^, 
and then the dynamics are very slow. 

One can also see that a strong magnetization is reached in the outer disc (/i ^ 1). This is very interesting 
because the condition for a powerful jet launching is satisfied at these radii. Then, the dynamic of the disc 
becomes considerably different (which is not taken into account in the code). If we consider the magnetic 
torque induced by the jet, we would assist to a new advection of the field from the outer edge of the disc. 
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Fig. 1. Density, accretion rate, magnetization and magnetic field in the disc at 3 different times. All the quantities are 
normalized. The solid line corresponds to initial state, the dash-doted line corresponds to 10 keplerian times at the outer 
radius, and the dashed line to 100 keplerian times at the outer radius. 



4 Conclusions 

The new advection of magnetic field supposed to occur at the end of the simulation would make a transition 
with the diffusion of the magnetic field observed in the simulation. Thus, we could assist to a phenomenon of 
magnetic tide. In this sense, this result is very encouraging for the study of transition from high soft to low hard 
state in X-ray binaries. However here, the accretion rate remains roughly constant during the simulation. But if 
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we work by imposing the accretion rate at the outer radius, a decreasing of it would enhance the magnetization 
at the outer radii, which is compatible with a transition between high soft to low hard state. 

Thank you to all the team from the laboratoire d'astrophysique de Marseille for the organisation of this conference. 
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